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Abstract 

For classical lattice systems, the Dobrushin-Lanford-Ruelle theory of boundary con- 
ditions states that the restriction of a global equilibrium state to a subsystem can be 
obtained as an integral over equilibrium states of the subsystem alone. The Hamil- 
tonians for the subsystem are obtained by fixing a configuration for the variables in 
the complement of the subsystem, or more generally, by evaluating the full interac- 
tion Hamiltonian with respect to a state for the complement. We provide examples 
showing that the quantum mechanical version of this statement is false. It fails 
even if the subsystem is classical, but embedded into a quantum environment. We 
suggest an alternative characterization of the local restrictions of global equilibrium 
states by inequalities involving only local data. 
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1. Introduction 



It is well-known that basic phenomena in statistical physics such as phase transitions, 
critical behaviour, and symmetry breaking cannot be modelled by finite systems as 
studied in classical or quantum mechanics. On the other hand, infinite systems can 
rarely be handled directly by finite computations. Therefore it is a fundamental task 
of statistical mechanics to develop techniques by which properties of infinite systems 
can be inferred from a study of finite subsystems. For the equilibrium statistical 
mechanics of classical systems on a lattice there is a standard tool for performing 
precisely such a step, namely the Dohrushin-Lanford-Ruelle (DLR) theory of bound- 
ary conditions [Isr] . It says that the effect of the ambient infinite system on a finite 
subsystem can be parametrized adequately by a choice of "boundary conditions" for 
the finite subsystem. Here, boundary conditions are identified with configurations 
of the particles outside the finite subsystem in consideration. They can then be 
used to classify the pure phases, or to select a particular pure phase in the event 
of a broken symmetry. In specific models such as, e.g., the two-dimensional Ising 
model, a small selection of boundary conditions suffices to obtain all phases. 

The notion of boundary condition as a configuration in the boundary region of a 
finite subsystem is the starting point both for proofs of existence of a phase transition 
and proofs of its opposite, the uniqueness of the eqiulibrium state. An example of 
the first is the Peierls argument [Sew] where one shows that the configuration of 
a small subsystem (a single spin), strongly depends on the configurations of the 
boundary region. An example of the other usage is the Dobrushin-Shlosman unicity 
argument [DS]. In this case, the influence of variations in boundary conditions is 
proven to become negligible when one surrounds a small system by larger and larger 
intermediate regions. 

It is one aim of the present paper to show, by way of explicit examples that, 
while the basic statement of the DLR-theory has a straightforward generalization 
to the case of quantum lattice systems, this generalization, along with some of its 
natural modifications, turns out to be false. But our intention is not merely to 
disprove a single statement. Rather, we believe that a careful study of examples 
may be helpful for the development of the much needed tools of quantum statistical 
mechanics which would be capable to serve some of the functions of the classical 
DLR-theory. The energy-entropy inequalities [FV] described in Section 4 are a 
promising step in this direction: they form a set of strictly local constraints on the 
restrictions of possible equilibrium states, which approaches an exact characteriza- 
tion of equilibrium states as the local region under consideration becomes larger and 
larger. For the restrictions to classical subsystems, an efficient evaluation of these 
conditions is possible (see Section 5), but, for the case of general quantum systems, 
better techniques of extracting information from these inequalities are still to be 
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developed. 



The term boundary condition, especially as used in the context of quantum 
statistical mechanics, carries several meanings, often not too clearly defined. "Free 
boundary conditions" refer to a special way of constructing finite volume Hamil- 
tonians from an interaction potential, namely as the sum over all potential terms 
localized in the given region. There are always many different interaction potentials 
representing the same infinite volume Hamiltonian ("equivalent potentials" [Isr]), 
which in this language lead to different "free boundary conditions" for the same 
Hamiltonian. Some properties, like the presence or absence of frustration [Mie, 
GW] , or the invariance of local Hamiltonians under the action of a symmetry group 
or quantum group [PS], depend very sensitively to this choice. "Periodic boundary 
conditions" arise when a lattice is shaped into a torus. The translation symmetry of 
the infinite lattice is then approximated by cyclic shifts. Another type of "boundary 
condition" is obtained by adding an overall small external field, which is removed 
afterwards [Ger] . Such fields are usually introduced for driving the system in a pure 
phase of broken symmetry. Typically, all these additional terms to the Hamiltonian 
leave the thermodynamic quantities unchanged. In the description of continuous 
systems "boundary conditions" are used to specify the Hamiltonian as a self-adjoint 
operator [BR]. This can be used to confine the system to a particular region, or 
to introduce hard-core repulsion. In principle, this is equivalent to adding strongly 
repulsive potentials to the Hamiltonian. Finally, conditions that select particular 
ground states in situations with degeneracy are also called boundary conditions. An 
example of such usage occurs in valence-bond-solid models [FNW] . 

For simplicity, we will only consider finite systems in this paper. The difficul- 
ties that we demonstrate for describing the local restrictions of equilibrium states 
will certainly not be lifted by considering infinite systems. (This could only hap- 
pen if more structure, like translation invariance, is incorporated into the scheme). 
However, the notions that we introduce for finite systems can be generalized in 
a straightforward manner to the infinite systems obtained in the thermodynamical 
limit. The adequate tools for doing so arc the DLR-equations for classical lattice sys- 
tems, the KMS-condition [BR] for quantum spin systems, and the energy-entropy 
inequalities for both. The basic input, namely the interaction between the lattice 
spins, now translates into a relative Hamiltonian for the classical systems, and a 
derivation that generates the dynamics for quantum spin systems. There is an, ap- 
parently local, description of equilibrium states for quantum spin systems known 
as the Gibbs condition [AI] (see also Sect. 6.2 of [BR]). Unfortunately, it involves 
a variation over all global extensions of the local state, and is therefore of little 
practical help. 

The paper is organized as follows: in Section 2 we briefly review the local 
"DLR-characterization" of equilibrium states of classical systems in terms of bound- 



3 



ary conditions. Section 3 introduces the corresponding notions for quantum spin 
systems. In Section 4 we recall the characterization of equilibrium states in terms 
of energy-entropy inequalities, while Section 5 examines the consequences of these 
inequalities for the restrictions of global equilibrium states to a classical subsystem. 
Section 6 contains the major part of our results, a succession of examples showing 
the failure of the direct quantum analogue of the DLR-theory. Where feasible we 
also discuss the consequences of the energy-entropy inequalities. 

2. Boundary conditions for classical 
equilibrium states 

In this section we will only consider finite, discrete classical systems, described by 
finite configuration spaces. We will concentrate on a distinguished subsystem, called 
the inner system, with configuration space X . The configuration space of the global 
system is X x y , where Y is the configuration space of the outer system. We assume 
that we have a complete knowledge of the contribution of the particles in the inner 
system to the total energy iiftot, meaning that we write 



The function {x,y) ^ X x Y ^-^ H{x,y) is the basic information that is at our 
disposal, whereas the size of the outer system and the interaction V is supposed to 
be less well-known. The minimal choice of outer system is fixed by requiring that 
H actually depends on the configurations in this minimal exterior region. More 
precisely we say that y and y' in Y are equivalent iff for all x E X, H[x, y) = H{x, y'). 
The quotient of Y by this equivalence relation is called the set of the boundary 
configurations Ibd and it is the minimal exterior configuration space that has to be 
considered. We will sometimes refer to it as the boundary system. 

The equilibrium state of the global system at inverse temperature (3 = 1/kT is 
given by the canonical Gibbs distribution, which assigns the probability 



to the configuration {x, y). The normalization factor is the canonical partition func- 
tion Z = V Q-l3Htot{x,y) ^ Unless we explicitly discuss temperature behaviour, we 
assume that /3 is included in the Hamiltonian, and therefore put /3 = 1. 

The restriction fXi^t of to X is easily computed. It assigns to x G X the 
probability of {{x,y)\y G Y}: 



Htot = H + V 



H fixed function on. X xY and 



(2.1) 



V only depending on Y 




(2.2) 






^-Htot{x,y) 



y€Y 
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where 



//J'(a;) = e-^'(^)/Z^ 



(2.4) 



is the Gibbs distribution on X defined by the averaged Hamiltonian 

Hy{x)=H{x,y) , 



(2.5) 



and where the X{y) are positive numbers adding up to 1. Thus formula (2.3) says 
that the global equilibrium state, reduced to the inner system X, is a convex com- 
bination of equilibrium states /i^ of the inner subsystem for averaged Hamiltonians 
Hy. The configurations y G Ybd can be viewed as boundary conditions, defining 
the averaged Hamiltonians Hy. Formula (2.3) solves the problem of giving a simple 
characterization of the equilibrium distributions, reduced to the inner system, in 
terms of the basic interaction H. It is a basic tool in proving both absence and 
occurrence of phase transitions. 

It is easily seen that all the Gibbs distributions are themselves limits of 
restrictions of equilibrium distributions of the global system to X . Indeed, it suffices 
to consider the minimal system X x Fbd and to take V — Xcy (where the function 
Cy takes the value 1 on y and zero elsewhere), and consider the limit A — > — oo. 
This makes the configuration y e Ibd infinitely attractive to the outer system, and 
yields a sequence of equilibrium states, whose restrictions converge to This 
procedure is reminiscent of the way one finds the Dirichlet boundary conditions for 
the Schrodinger equation in a finite region, by introducing larger and larger confining 
potentials. 

The set of restrictions of global equilibrium states /i, determined by Hamilto- 
nians Htot satisfying (2.1), is convex. This is easily seen as follows: if Y and Y' are 
two outer configuration spaces for the subsystem, consider then Y" — Y UY' and 
the total Hamiltonian H^'^^ = H + V'\ where V"{y) — V{y) if y belongs to Y and 
V"{y) = V'{y) + C else. By varying the constant C e H, we can produce any con- 
vex combination of the restrictions of the equilibrium distributions corresponding 
to -fftot and -fftof Taking the conclusions of the last three paragraphs together, we 
arrive at the following statement, whose quantum versions form the subject of this 
paper: 

for a measure ^ on X the following two statements are equivalent: 

(1) fj, can be approximated by the restriction of Gibbs distributions for 
Hamiltonians in the class (2.1). 

(2) jJL is a convex combination of the Gibbs distributions , formed with 
Hamiltonians Ry (2.5), averaged with one outside configuration y e 
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Especially suggestive in this equivalence is the identification of the "extreme" el- 
ements of the set so described: in terms of the reduced equilibrium distribu- 
tions (1), these are obtained by making the outside potential V very large, whereas 
in (2) the same elements play the role of the extreme boundary in the sense of convex 
sets. 

This suggestive picture loses a bit of its appeal when we consider more general, 
less "extreme" averages of the Hamiltonian H: rather than evaluating H on a single 
external configuration y E H, we can also consider 

HP{x)= ^ p{y)H{x,y) , (2.6) 

where p is a probability distribution on Ybd- Such boundary conditions arise nat- 
urally when considering random systems. They will also appear when we couple 
classical systems to quantum outer systems. One can easily find examples in which 
the equilibrium states of such averaged Hamiltonians are not contained in the 
convex set described above (see Example 1 for an illustration). 

Besides equilibrium distributions, ground states can be considered as well. A 
distribution is a global ground state if it minimizes the energy functional ji i— > 
li{H^oi). As this functional is obviously affine, the ground states form a face of 
the simplex of probability measures oiv X xY , and the same is of course true for 
their restriction to the inner system X. By an argument quite similar to that for 
equilibrium states, one verifies that the set of reduced ground states, when Y and 
V vary as in (2.1), is still a face of the probability measures on X. A configuration 
xq & X is an extreme reduced ground state iff there is a j/o £ ^bd such that 

H{x, yo) > H{xo, yo) for all x e X . 

The situation is quite analogous to that for equilibrium states: in order to char- 
acterize reduced ground states, it suffices to minimize the energy of the averaged 
Hamiltonians (2.5) where y runs over the configurations of the boundary region 
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3. Extension to quantum systems 



There are two essential properties of classical systems that underlie the description of 
equilibrium states of global systems in terms of local Gibbs distributions. Firstly, all 
probability distributions on a composite system are convex combinations of product 
measures, namely the point measures described by pairs of configurations. This 
allows conditioning as in equation (2.3). Secondly, for functions / and g, exp{ f+g) = 
exp / exp g. Both properties totally fail when we pass to quantum systems. 

Quantum systems are described in terms of algebras of observables rather than 
configurations. The observables of a discrete, fully quantum mechanical system form 
a matrix algebra, say, the algebra Aid of complex (i x (i-matrices. On the other hand, 
the observables of a discrete classical system with finite configuration space X are 
the complex- valued functions C{X) on X. This is an Abelian algebra, which can be 
identified with the diagonal dxd matrices, d being the number of points in X. In this 
paper, we only consider finite-dimensional algebras of observables. All such algebras 
decompose into sums of matrix algebras and may also describe quantum-classical 
hybrid systems. Probability distributions are expressed in algebraic language as 
linear functionals, called states, which assign to each observable its expectation 
value. On a finite-dimensional algebra A all such functionals are of the form 



where the density matrix D,^ is positive, and TrD;^ = 1. 

We consider systems split into an "inner" and an "outer" subsystem. On the 
level of the algebras, this translates into a tensor product structure of the global 
algebra of observables A — Aint ® ^ext^ the elements of A can be written as finite 
linear combinations of elementary tensors A<Si B, A e A, B e B. A t-^ A ® B and 
B ^ A (S) B are linear and the multiplication in A is given by 



As in the classical case, we will assume that we have a complete knowledge of the 
interaction of the particles in the inner system with the exterior. More precisely 



Here, ^bd denotes the smallest *-subalgebra of Aext such that H e ^int ^bd- 
For classical systems ^bd is precisely C(Ybd)- The equilibrium state u of the global 
system is given by the Gibbs density matrix 



u;{A) = Tr{D^A) 



{Ai ® Bi) {A2 ® B2) =^1^2® Bi B2 



H e ^int <8) .4bd is a fixed Hermitian operator and 
H,^t = H + l®V, V = V* eAext^A^A . 



(3.1) 




/2 



(3.2) 
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The normalization factor Z is the canonical partition function. We can restrict the 
state LO to the inner system to obtain the state los^a on .Aint: 



a;int(^) =a;(^® I), ^ G Ant • (3.3) 

The density matrix of a;int is the partial trace of the one for lo with respect to Aext • 
We now define the set of reduced equilibrium states 



= I'^int ^ equilibrium state for ii^tot ^ in (3-1) | ) (3-4) 

where the bar denotes the closure. The set 5r is completely determined by B. and 
we will sometimes explicitly refer to that dependence by writing = Sji{H). By 
the same argument as in Section 2, £r is a convex set and we are interested in the 
following question: can one give a simple description of £r in terms of HI 

There are two natural ways of averaging the Hamiltonian H: either against 
pure or against arbitrary states p on Ahd (compare with (2.5) and (2.6)). In either 
case, the averaged Hamiltonian is obtained as 

W{H) = Y,p{B,)A, , 

i 

where H = Ai ® Bi with Ai G Aint and Bi G Ahd C Aext ■ The sets of Gibbs 
states on Aint corresponding to these averaged Hamiltonians will be denoted by 

£P^^«= = jo; I equilibrium state for E^(iy) , p pure | (3.5) 

and 

^AH = equilibrium state for JEF{H) , p arbitraryj . (3.6) 



In [Isr] , the set of pure states on Ahd is considered as a candidate for quantum 
boundary conditions. As in the classical case, this choice can be made plausible by 
considering perturbations of the Hamiltonian by large potentials: let P be a minimal 
projection in Aext and consider the Hamiltonians i?tot('^) = H + Xl<Si P- Note that 
P, taken as a density matrix, also determines a pure state p on Ae^t- We claim that, 
as A ^ — oo, the restrictions of the equilibrium states for ii'tot(A) to Aint converge 
to the equilibrium state for lEf (i/). In order to verify this, we consider \~^H as a 
small perturbation of — AI(8)P, obtaining eigenvectors il>ai^)i and eigenvalues rja{\), 
which depend analytically on A~^. Note that the eigenspaces of I ® -P are highly 
degenerate, and we must chose the basis in 0**^ order such that ?/'a(oo), a = 1, . . . , (i 
is the eigenbasis of the operator (I (g) P)H{1 P) and the remaining iIjo,{oo) are an 
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eigenbasis of (1 (I — P))H{1® (l — P)). Let r]'^{X) denote the eigenvalues of these 
operators. Then the equihbrium state a;'*' is given by 



The exponent is 

Ar7a(A) 



-X + f]'^{X) + 0{X-^) fora = l,...,d 
77^(A) + 0(A-i) fora>d 



Hence the terms with a > d are neghgible in the sum and, after canceUing the factor 
between the sum and 2{X), we can take the hmit A — > — oo, obtaining 

1 

u^{A) = —— V e-<(-)(V'«(oo), ^V«(oo)) 

= (l/Z)Tre-(^®^)-^(^®^)yl . 

An infinite system generalization of this result can be found in [RW]. Now, if P 
is the one-dimensional projection onto the vector the 1-eigenspace of (I ® P) 
consists precisely of the vectors ip ® (p where ip is in the Hilbert space of the inner 
system. By identifying ip0ip with ip, the operator (I(8)P)i7(I®P) becomes M'{H), 
hence u>°° is the Gibbs state with this Hamiltonian, as claimed. 



As any state on can be extended to a pure state on ^ext for some suitable 
embedding of into a ^ext; it follows that 

couv^ahC^r . (3.7) 

Obviously, C £au and the states in can be obtained by restricting in the 

construction of above to the case ^ext = «4.bd- 



In the classical case we found that 8^ = conv Sp^'^. In Section 6 we will examine 
whether the converse inclusion of (3.7) holds, i.e., whether we have 

f^RCconvf^AH • (3.8) 

We will refer to this as the DLR-inclusion. By (3.7) it is in fact equivalent to the 
equality £^r = couv^'ah- 



A state u) of the global system is a ground state of i^tot if it minimizes the energy 
uj{Htot)- For a fixed -fftot, the ground states form a closed convex set, which is even 
a face of the state space. Restricting a pure state of the global system ^int ^ext 
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to the inner system will in general destroy its purity. The set ^r, of reduced ground 
states is similar to £r 



uj ground state for i^tot as in (3.1) | . (3.9) 

It is a closed convex set that contains the convex hulls of the sets ^Xh^® ^ah 
which are defined as the sets of ground states on Aint with respect to the averaged 
Hamiltonians IE''(iy), where p varies either over the pure or over the general states 
on ^bd- 




4. Energy-entropy inequalities 



An equilibrium state can be characterized as a solution of the variational principle 
of thermodynamics: 

uj ^ F{uj) = U{uj) - S{uj) 

attains its minimum at the Gibbs state. F is the free energy, it is the difference 
between the internal energy U{u3) and the entropy S{uj) = —Tt D^jlnD^^. It is 
possible to characterize equilibrium states by a set of inequalities that express in a 
differential form that the Gibbs state minimizes the free energy [FV]. Consider a 
quantum system described by the dx d matrices Aid and let H be the Hamiltonian 
of the system. For any observable A e Md, the Gibbs state oj satisfies the "energy- 
entropy inequalities" 

u;{A*[H,A])>u;{A*A)ln^^^^^ . (4.1) 

Conversely, if (4.1) holds for any choice of A G Md, then a; is the Gibbs state defined 
by H. The main advantage of these inequalities, which express the energy-entropy 
balance in equilibrium, is that they involve the Hamiltonian in a linear way. 



Consider now again the situation of Section 3. Any equilibrium state iv of the 
global system will satisfy the energy-entropy inequalities (4.1). In particular, the 
inequalities must hold for A G Aint and therefore: 

u;{A*[H,,„ A]) = uiA*[H, A]) > uiA*A) In , ^ G Ant • (4.2) 

Note that the left hand side involves only the restriction of a; to ^int®«4.bd5 whereas 
the right hand side involves only the restriction to .Aint- In either case, this condition 
is "local" in the sense that it does not require knowledge of the global state on ^ext- 
Moreover, if this condition is satisfied for ^int varying over the observable algebras 
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of arbitrary bounded regions, the inequalities (4.1) hold on the whole algebra, and 
hence u> must be an equilibrium state. Thus at least some of the key requirements 
for a useful replacement of the DLR-equations in the quantum case are satisfied by 

We denote by See the set of states on ^int defined by 

^EE = |<^| 3 a;~ extending uj and satisfying (4.2) | 

The set S-^e is a closed convex set and it contains Indeed, in order to show 
convexity, let u and u' be states in See- We can then find extensions lj~ a-nd u;'~on 
^int ® Ahd that satisfy 

a;W~(A*[if,A])>a;W(A*A)ln^^|^, A G Ant • 

For A e [0, 1], Xuj~+ (1 - A) extends Aa; + (1 - A) u;' and for A e Ant: 

(Aa;~+ {1-X)lj'~){A* [H, A]) 

= Xco~iA* [H,A]) + {1-X) to'- {A* [H, A]) 

= - h.(A'A) m 4^ + (1 - A) JiA-A) m 44^) 

V ^ ' uj{A*A) ^ I \ I uj'{A*A)) 

Xu{AA*) + {l-X) oj'{AA*) 



> -{Xuj{A*A) + (1 - X)J{A*A)) In 
= {Xuj{A*A) + (1 - A) u\A*A)) In 



Aa;(AM) + (l-A) uj'{A*A) 
Xuj{A*A) + {l-X) u\A*A) 



Xuj{AA*) + {l-X)u'{AA*) 



There are versions of the energy-entropy inequalities for classical and hybrid 
systems too. The inequalities (4.1) are of no use in this case because in a classical 
system all commutators vanish and A* A = AA*. Hence (4.1) is trivially satisfied 
as > 0. We therefore need a more general type of inequality that also covers 
hybrid systems. Consider, e.g., a coupled classical-quantum system described by 
the observable algebra C{X) ® A. If G X denotes the function which is 1 at a; 
and otherwise, we can expand the observables A of this system as 

A=^ea;® A{x) , (4.3) 

x€X 
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where the A{x) are elements of A. Thus observables are identified with ^-valued 
functions on X, and the algebraic operations are defined pointwise in x. In par- 
ticular, A is Hermitian iff A{x)* = A{x) for all x, and the identity element is the 
function which is equal to I e ^ for all x. A state oj on such a system is of the form 

a; = ^ (5a; (g) a;^ , (4.4) 
xex 

where 5x is the evaluation at x and the are positive functionals on A. The 
expectation of an observable (4.3) in the state (4.4) is computed as 

oj{A)=J2^^iA{x)) . 
With the identity element for A, we find the normalization condition 



The useful energy-entropy inequalities in this situation are 

u;-[A*{H{y)A-AH{x)))>u;^{A*A)ln^^^^^^ , (4.5) 

x,y & X and A & A. Prom this, we find restricted energy-entropy inequalities 
analogous to (4.2). Assuming that the inner system is the classical subsystem, i.e., 
Aint = C{X), we have to evaluate (4.5) for A = al. Thus AA* = A* A, and a factor 
I a I can be cancelled from (4.5). Hence, for a classical inner system, the suitable 
definition for See is the set of all probability measures n on X for which there exist 
non-negative functionals u>^ on A with ij,{x) = u>^{T) and such that 

u;^{H{y)-H{x))>^{x)ln^, x, y e X . (4.6) 



12 



5. Systems with a classical interior algebra 



With the exception of Example 3, aU examples in Section 6 will be "hybrid systems" 
with a classical inside system with n = finitely many configurations, i.e., we have 
«4.int = . . . , n}). The interaction is given by a Hermitian element H e Aint ® 
vAext, i-e., according to (4.3), by a Hcrmitian-valued function x H{x) G ^ext- 
The subalgebra C ^ext is the algebra generated by I, and all H{x). The 

perturbation l^V by a potential V — V* & Aext is represented by the corresponding 
a;-independent function. The reduced equilibrium states of the total Hamiltonian 
i^tot = H + l<SiV are then the probability distributions 

f^U^) = ^T^{e-(H{x) + V)^j ^=i,...,n , (5.1) 

where Z is the normalization factor making this a probability distribution. If p 
is a state on ^ext; the averaged Hamiltonian 'Wf{H) is the function h = lE^(iJ) e 
C({1, . . . , n}) taking the value h{x) = p{H{x)) at x. (Note that only the restriction 
pf^bd enters this expression). This leads to a probability distribution of the form 

^W{^) = ^e-P^H{x)) ^ = 1,...,^ . (5.2) 

General states ui on Aint ® Aey± are given by (not normalized) positive functionals 
uj^ on ^ext) as in (4.4). By (4.6), the energy-entropy inequalities for such a state 
mean that 

u-(H{y)-H{x))>p^^{x)\n^^^^ , a;,y=l,...,n . (5.3) 

^ ^ /iEE(?/) 

Of course, the restricted state /xee G £ee is related to the functionals uj^ by 
A*EE(a;) = <^^(I), with uj satisfying this inequality. 



It is sometimes convenient to work with the logarithms of probabilities, i.e., to 
describe a state p on ^int by the n — 1 numbers 

(ln//)(x) = ln(//(x)) - ln(//(l)) , for x = 2, . . . , n . (5.4) 

The image of the various sets £^xx C 1R"^~^ under this map will be denoted by 
In^xx- In this logarithmic scale Sab. becomes 

hi^AH = {^elR"-^|^, = p(if(l)-ff(a;))} , (5.5) 

where p runs over all states of ^bd- This is an affine image of the state space of 
vAbd- Note that this set is always convex but, due to the non-linearity of exp, £ah 
itself is not (compare Figure 1). 
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When ^int is classical, we can also give a fairly explicit description of See- Since 
we want to interpret (5.3) as a condition for the probability distribution ij,ee{x) = 
c<;^(I), it is convenient to write the positive linear functionals a;^ determining u> as 
u>^ = heb{x) ■ p^, where the are now states on ^ext- Note that no HEE{y) can 
vanish: since u; 7^ 0, some Hee{x) must be non-zero, and by (5.3) ln{pEB{x) / lJ>EE{y)) 
is not infinite. Hence we may divide (5.3) by Pee{x) and we find 

lnpEEiy)+p''iH{y)) > \iip.EE{x)+p%H{x)) . (5.6) 

In these terms, the energy-entropy condition on //ee is that, for every x, there 
is a state on ^ext such that (5.6) holds for all y. We consider first the case 
X = 1. Then, using the equations (5.4) and (5.5), we conclude that there is a vector 
^ G In^AH, namely = p^{H{l) — H{y)), such that In //ee — C li^s in the cone Ci of 
coordinatewise positive vectors in WC^~^. In shorthand notation, this is written as 
In jjEE G ln£^AH + C11 where, as usual, the sum on the right is to be interpreted as 
1^ + ^ C £ In^AH, 7/ e Ci|. The conditions for a; > 1 can be written in a similar 

form after subtracting ln//EE(l) + P^{H{1)) from both sides of the inequality when 
X ^ y. Thus, with the cones 

^ ^ . UelR--^|Vy^,>0} x^l 
{^elR"-^|Vyey>exand^,<0} x>l 



we get the formula 



lnSEE= OilnSAK + C^) . (5.8) 



x=l 



This construction is illustrated below in Example 5 (see Figure 4). An amazing 
feature of this set is that both In See and its exponentiated version See &re convex 
(see Section 4 for the second statement). This seems to be in confiict with the 
non-linearity of the map In . However, this conflict is resolved by considering the 
special form of the cones (5.7): the cone Cx corresponds to that subset of the state 
space in which \np{x) < \np{y) for all y. Such inequalities survive exponentiation, 
i.e., the cone Cn^ +\n n at n corresponds in the state space to the set of probability 
distributions p' such that p'{x)/p{x) < p'{y)/p{y) for all y, which is a convex set. 
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6. Examples 



Since our main aim is to demonstrate the failure of the DLR-inclusion in the quan- 
tum case, our main results are in the form of examples demonstrating the success 
or failure of these ideas in various situations. We have chosen a succession of ex- 
amples to illustrate various degrees of failure of the DLR-inclusion £^ C couv^^ah- 
Mathematically, the strongest result is the last example, in which both the inside 
system .Aint, and the boundary algebra Ahd sixe classical, and Aext is the algebra 
of 3 X 3-matrices. However, this example is rather indirect, so we will give simpler, 
direct examples as well, hoping to help the reader to develop an intuition of what 
exactly goes wrong. We emphasize that, while we treat only finite systems, our 
results are equally valid for infinite outside systems. In fact, every example that we 
give can easily be enlarged to an example with infinite outside. In the same spirit 
we are looking mostly for examples with classical A-mt- usually these can be ex- 
panded easily to examples with non-commutative Aint- Moreover, they are simpler 
to treat and, being closer to the classical DLR-situation, they bring the failure of 
the quantum generalization into sharper focus. In producing the examples we found 
the help of a computer algebra program [Mat] a very valuable tool. 

The symbols in the heading of each example are a shorthand for the type of 
algebras used. Thus "Q5C3" is an example with a quantum inside with Aint — -Mb, 
and an interaction H contained in a subalgebra of Aint ® -^bd, which is isomorphic 
to Aint 2, 3}), i.e., the inside system interacts with 3 classical configurations. 

Of course, the whole algebra Aext is never assumed to be classical, since then the 
DLR-theory would simply be valid. If the number of configurations or Hilbert 
space dimensions is irrelevant, wc just write "n". Each example begins with a short 
description of the model and of the claims about various inclusions to be seen in 
this example (set in italics), followed by the verification of these claims and some 
additional remarks. 



Example 1: C4C3; ^ext classical 

With this example, where Aext is Abelian, we demonstrate that cony Sab. 't- conv^^^^^ 

It also serves to illustrate the lack of any definite convexity properties of the map 
which takes Hamiltonians to Gibbs distributions (see Figure 1). Explicitly, the in- 
teraction is given by the matrix H{x, y), x — 1,2,3, 4, y = 1, 2, 3, as 



The numerical evaluation in this example is straightforward, and the result is shown 



H = 
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in Figure 1. We include this example mainly to illustrate a remark made in Section 2, 
namely that, due to the non-linearity of the exponentiation map, general averaged 
Hamiltonians (2.6) do not give equilibrium states inside the DLR-triangle £^r, which 
is also marked in the figure. 



Example 2: C2Qn 

When the inside system is a single Ising spin and the outside is arbitrary, the DLR- 
inclusion holds, namely 



The Hamiltonian is specified by two Hermitian H{1),H{2) G Ain- Then, by (5.5), 
ln£^AH is just the set of numbers p{H{l) — H{2)), where p ranges over all states 
of ^ext- Thus In^AH is the interval [r7_,?7_|_], with endpoints given by the largest 
and smallest eigenvalues of H{1) — H{2), i.e., the smallest r]+ and the largest rj- 
such that r]-l < H{1) — H{2) < ry+I. Moreover, the same interval is obtained if we 
restrict the states p to be pure, i.e., Sau = ^ah'' ~ [V-jV+]- 

On the logarithmic scale, p^^^ e £r from (3.4), is given by the single number 

r]{V) = InTr exp -{H{2) + F) - InTr exp -(iy(l) + V) . 

Now the function A i-^ InTr exp(— A) is monotonely decreasing with respect to 
the operator ordering for A. Moreover, InTr exp(— A — 77I) = InTr exp(— A) — ry. 
Since H{2) < H{1) — r/_I, we get rj^V) > rj-. Similarly, we get the upper bound 
r]{V) < ?7-|-. Hence ^r, C [r]-,rj+] = S^^, and, from the general arguments in 
Section 3 concerning large potentials, we get the inclusion ^ah'' hence 
equality. 



Example 3: Q2Q2 

We now replace the Ising spin of the inside system by a quantum mechanical spin-^, 
i.e., we take Aint = M.2- For the outside, we will also have Aext = -^.2, and define 
the Hamiltonian by 

i?(Q;, 7) = (Ji (g) (Ji + (72 (g) (72 + 7 0-3 (g) 0-3 - a 0-3 (g) I , (6.1) 

where a > and 7 > are parameters, and Uj, z = 1, 2, 3 are the usual Pauli matri- 
ces. We claim that, in this model, the reduced equilibrium state tof^^ e S^_(^PH{a, 7)), 
where the superscript denotes that the potential V is vanishing, is not contained 
in conv £Aii{PH{a,'y)) . Hence the DLR-inclusion fails in the following cases: 
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(1) /3 = 3, a = .3, 7 = .05 (see Figure 2) 

(2) 7 = 0, and all 13 > Q, 

(3) for < 7 < .8, and j3 = oo, i.e., for ground states. 

What makes this operator tractable is that it commutes with SU(2)-rotations around 
the 3-axis. In the joint eigenbasis of (J3 I and I (8) cts it corresponds to the matrix 



if (a, 7) = 



/-a + 7 0\ 

-a -7 2 

2 a-7 

V ^+7/ 



The partition function for if (a, 7) and the expectation of (73 (8) I are readily com- 
puted: 

Z(/3,Q!,7) = Tre-^^^^^'T') = 26"^^ cosh(/3a) + 2e^^ cosh(/3i?i) 
1 9 

^int(^3) ^^-Q^^"^ ^(/^' ^) 

^g/3(a+27) /g2/3i?i _ l^ ^ ^ /3i?i /g2a/3 _ ^\ 

^ ' ' ^ ' (6.2) 



i?l (e/5(«+27) + e^^i + e'3(2a+^^l) _^ g/3(a+27+2i?i)) 

where i?i = vTTc?. Since H{a,'^) commutes with rotations around the 3-axis, so 
does a;piit, and hence it is completely determined by the expectation of (73. 

If p = I ( 1 , With e R and a;^ < 1, is a 2 x 2-density matrix, 

we find 

]E^(ii(a,7))=f-" + :^^3 ^^-^^2 

Y xi + ia;2 OL — 7x3 

The equilibrium state a;^jj of lE''(ii(Q!, 7)) can be evaluated with the formula 

e-/3^ 1^ tanh/jy/d^ ^ 



Tie-(^H 2 2x/d^ 

for ii a traceless 2 x 2-matrix. £^ah is now a body in the state space of M2, 
parametrized by the triple (xi, X2, X3) determining p. It is the image of an ellipsoid 
in the space of 2 x 2 Hamiltonians under the non-linear map (6.3). This map 
preserves the symmetry under rotations around the 3-axis, hence £^ah will also have 
this symmetry. Figure 2 shows a section of Ear in the (l,3)-plane of the state 
space. It has to be completed by rotating it around the vertical axis, thus making 
£^AH the shape of a mushroom. The point marked on the axis is o'Pnt- clearly 
not contained in the convex hull of £r,. 



17 



At high temperature (small P) one can expand all exponentials, and one finds 
that, in first order in (3 and as long as 7 ^ 0, the convex hull of Sau does contain 
a;Pjt. However, if 7 = 0, the first orders coincide, and we get 

<uM - <t(^3) > 3 (2 - ^1 - ^2) + 0{(5') > . 

Hence the DLR-inclusion fails even for small /5, and it can be seen numerically that, 
with 7 = 0, it fails for aU /3 > 0. 

In the opposite direction, as (3 grows, the points in £ah move radially out to the 
extreme boundary of the state space and, in the limit, we obtain the pure ground 
states of the Hamiltonians Wf{H{a,^)). Hence, for /5 = 00 



^Ah(^3) = / . o , > 



where the second expression is the minimum of the first with respect to variation 
of the Xi, attained for x\ + = — x'^ and x-i = 7/0;. The limit /? — > 00 of a;-^^ 
can be obtained directly from (6.2): of all the exponents appearing in this equation, 
/?(« + 27 + 2i?i) is clearly the largest, so the limit is a/Ri, or 

Hence, for all p, we have 

'^ah('^3) > ^intl^a) for /3 = 00 and < 7 < . 

Thus the DLR-inclusion fails in general for quantum ground states. For this conclu- 
sion it is crucial that both systems are non-classical. In our further examples we will 
always take ^int Abelian. The pure states of the tensor product ^int ®-^ext, such as 
the extremal ground states of H, are then product states and can be found by min- 
imizing W{H) for each p. Thus, for classical Ant) the ground state DLR-inclusion 
holds. 



Example 4: QnCm; = >^bd ® ^ 

The DLR-inclusion S]^ = convf^^'^^'' holds when the boundary algebra is purely clas- 
sical and when the exterior algebra Aext factorizes into a tensor product C{Y\)(^)®B, 
B arbitrary. 

From the general argument of Section 2 we always have conv S'^^ C Conversely, 
with Htot = H -\- 1® V, H and V functions from Fbd to ^int and B, we have for 
any A e Ant 

E TiA..«B(e-<^ + ^)^)= E li«(e-''(^))^A..(e-^ffe)^) . 
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But this precisely means that Sji C convS^^. From Example 1 it becomes clear 
that generally conv^^'^'^'^ is strictly contained in conv^AH- Therefore the inclusion 
conv£^AH C £r, that holds if we allow for general embeddings of C(Ybd) in an exterior 
algebra ^ext, is reversed if we restrict ourselves to product systems Ant®C(ybd)®'S- 

A special case of this example is the case where all of ^ext is classical. The rep- 
resentation of equilibrium states by conditioning with respect to a classical outside 
is reminiscent of the work [FF] . In that project one also considers quantum systems 
split into inner and outer region. The aim is then to represent general quantum 
states as integrals of states conditioned with respect to a classical subalgebra in the 
outside region. In the present, finite-dimensional setting, this is a trivial operation. 



Example 5: C3Q2; ^ext = M2 

By Example 2, the smallest non-trivial case of a classical inside algebra is Ai^t = 
2, 3}) and we stay with Aext = -M.^ i-e., we consider only the subset ^r(A^2) C 
Sji of reduced equilibrium states coming from potentials V e Al2- Thus H is given 
by three 2 x 2-matrices, which we choose as 

H{1) = , H{2)^^a3 , and i?(3) = ^cJi . 

In this example, we will determine the three sets Sah C £^r(A12) C £ee explicitly 
and show that both inclusions are strict (see Figure 3). 



The averaged Hamiltonians 'EF{H) depend only on the two expectations p{a'z) 
and ^(0"!), and therefore form a two-dimensional disc in the space of Hamiltonians. 
Note that, because p(a"2) does not enter, we get the same disc when Hamiltonians 
are averaged only by pure states p. The disc of Hamiltonians is (up to a factor) 
identical with ln£^AH (see the circle in Figure 3). 

In the logarithmic scale jJL^-^ is given by a point in the plane with coordinates 

r,(v\-h,uy - AnT¥exp(-if(2)-V)-lnTrexp(-V)\ 
VKV ) -^^l^K^- y^^rj^ exp(-ii'(3) - y) - InTr exp(-F) ) ' 

By choosing some random V , it is easy to convince oneself that in this example 
^k{,-M.2) 't- ^AH, i-e., we get points outside the circle in Figure 3. The precise 
determination of S^{M.'2) is more difficult. Consider the map V ^ ?7(^)- We know 
that, as A — > cxo, r\{XV^ approaches a point in £^ah- Therefore the boundary points 
of £^r(A^2) must be of the form ?7(F) with finite V , and at such points the map r\ 
must be singular, more precisely, the rank of the Jacobian must be one. It turns 
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out that this can happen only in the plane V = xai + zas, along three disjoint 
arcs. Each of these arcs extends to infinity, and describes a piece of the boundary 
of S-R_{J^2) connecting with the boundary of £ak- Solving for the critical points in 
the {x, 2;)-plane can only be done numerically and the result is shown in Figure 3. 

In this example we can also evaluate the energy-entropy inequalities (5.3) . Using 
the general method outlined in Section 5, we obtain Iu^'ee as the intersection of 
the three convex sets shown in Figure 4. The result is also included in Figure 3 for 
easier comparison. 



Example 6: C3C2 

In the previous example the algebra Ahd was non-Ahelian. We now move still closer 
to the classical DLR-situation by requiring this algebra to be Abelian too. We will 
retain the choice Aint =C({1,2,3}) of the previous example, but set Ahd =^({1,2}). 
More generally, the considerations of this example are valid if Amt — ^({1;2,3}) 
and if the three operators H{x) e .4extj determining the interaction, are of the form 

H{x) = hQ-\- axhi + bxl , 

for some Hermitian ho, hi G ^ext o^nd real constants a^, b^. We then find that the 
DLR-inclusion holds, i.e., = couv^ah- Moreover, the inclusion £r C £ee is 
strict (see Figure 5). 

Without loss of generality, we can take a2 > oi > a^. Then we can find a A, 
< A < 1, such that 

H{l)^XH{2) + {l-X)H{2,) + yl , (6.4) 

where y = bi — A62 — (1 — A)63. Hence, we have 

In //(I) = A In fx{2) + (1 - A) In ^(3) -y , for e £au ■ (6.5) 

Thus Sah hes on a line in the state space, which appears straight on the logarithmic 
scale. The set Sah is parametrized by the expectation p{hi) in arbitrary states p 
of ^ext and the extremal eigenvalues r]± of hi (i.e., the best constants with rj-l < 
hi < 77+ 1) determine the Sah as a segment of that line. An alternative description 
of the endpoints is by the equation of the line passing through them, or equivalently, 
by the coefficients C2, C3 of the linear inequality 



/i(l) > C2//(2) + C3//(3) 
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(6.6) 



which becomes an equahty precisely at these endpoints. Inserting the known end- 
points /i± and solving the linear system, we find 

g?7_(a2-oi)+r7+(a2-ai)-|-(62-6i) /'g?7_(a3-ai) _ ^r]+{a3-ai)\ 

Co — — 

g77+(a2-ai)+r7_(a3-ai) _ grj- (a2-ai)+r7+(o3-oi) 

gr7_(o3-oi)+r7+(a3-ai)-|-(63-6i) /'gr7+(o2-oi) _ ^r}-{a2-ai)\ 
g77+(a2-ai)+r7_(a3-ai) _ gr7_(a2-ai)+r7+(o3-oi) 

The only information we will need from these formulas is that, since ?7_|_ > r]- 
and 02 > oi > as, both constants are positive. This implies that, for all 77 with 
V- <V< V+, 

Indeed, this is an equality for rj — rj^ and rj — rj- and follows for the intermediate 
values by multiplying both sides with e"'^'^ and invoking the convexity of the expo- 
nential function. Since p{hi) G [?7_,?7+] for all states p on ^ext, we find that Sah 
is the set of probability distributions satisfying equation (6.5) and inequality (6.6). 
(See the thick line in Figure 5). 

On the other hand, since the interval [rj- , 77+] contains the spectrum of hi and 
both {H{1) - H{2)) and {H{1) - iy(3)) are functions of hi, we obtain from (6.4) 
the operator inequality 

I>C2e^W-^(2)^g^g//(l)-H(3) _ 

Combining this with the Golden-Thompson inequality [Ara], i.e., Trexp(^-I-S) < 
Tr(exp(yl) exp(5)), we get 

C2TVe-^(2)-i/^^_^^^g-//(3)-v 

= c2Tre(-^(^)"-^(^))"(-^(^)+^) +C3Tre(-^(^)"^(^))"^'^(^)+^) 

< Tr (C2 e^(l)-'f^(2) ^ gg gH(l)-H(3)^ ^-H{i)-v 

<TVe-^«-^ . 

Hence, all probability distributions /i e £r also satisfy inequality (6.6). On the 
other hand, we may use the convexity of the function 1-^ InTr exp(— yl), together 
with (6.4) to obtain that, for jj, e £r. 

In < A In ijl{2) + (1 - A) In //(3) - y . 

Thus we have shown the DLR-inclusion £y{, C Exr- The determination of £ee uses 
the method described in Section 5 and illustrated in Figure 4. 
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Example 7: CnC3 



The previous example might have nourished the hope that, at least for purely clas- 
sical interactions (Ahd Ahelian), the DLR-inclusion might survive. What we will 
show now, however, is that S^, <f. conv£^AH in this case. The following Proposi- 
tion reduces our question to a well-known failed conjecture, which was disproved by 
Gaudin [Gau]. 



Proposition. Let ^ext — -M.^, ^nd denote by the subalgebra of diagonal ma- 
trices. Suppose that, for every n and every Hamiltonian H e (S)T)s, 
the DLR-inclusion C conv Eau holds. Then, for any Hermitian 3 x 3-matrix V, 
there is a positive measure jiy on such that 

Trexp ^- ^0 o j - ]/ j = J /ly {dXi dX2 dXs) e~^^»^» . (6.7) 

Proof : Let us fix some Hermitian 3 x 3-matrix V. By Ei G I^s we denote the 
matrix with "1" in the i^^ place of the diagonal, and all other entries zero. For 
X e we set H{x) = Yl^=i^i^i- Consider a finite subset X C R^. Then H{x) 
defines an interaction H G C{X) (g) V^. The left hand side of equation (6.7) is 
/^reI-^)' X E X. On the other hand, the average of H with respect to a state 
on (or, equivalently, a state on V^) is characterized by three numbers Aj > 
with Aj = 1, and TEi^{H){x) — ^y the inclusion C Sau we can 

find a measure Hy with finite support, such that equation (6.7) holds for all a; G X 
with jjy for fiy. Now consider the net of finite subsets of X C IR"^, ordered by 
inclusion. We may assume that each X contains the origin, so that each one of the 
associated measures Hy is normalized to the same constant Tr exp{—V). Then, by 
weak*-compactness, we find an accumulation point fxy of this net of measures and 
it is clear that equation (6.7) holds for this measure and all a; G H. 

■ 

We may paraphrase this by saying that, given the validity of the inclusion 
C ^^AH, the expression Tr exp A, considered as a function of the diagonal matrix 
elements of A, is the Laplace transform of a positive measure. This is precisely the 
statement investigated by Gaudin [Gau] and proven to be wrong by analyzing the 
measure fxy, which is, of course, uniquely determined by V due to uniqueness of 
Laplace transforms. 

Rather than considering the expression Tr exp ^ as a function of its n diagonal 
elements, one can also study this expression along a single straight line in the space 
of Hermitian n x n-matrices. The following conjecture about this situation was 
formulated by Bessis, Moussa, and Villani in 1975 [BMV]: 
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Conjecture. Let A and B be Hermitian n x n-matrices. Then there is a positive 
measure jjL such that 



A, possibly signed, measure satisfying this equation always exists and is 
uniquely determined by A and B. Hence, the issue is only the positivity of this 
measure. Despite many attempts, this conjecture is still open, even for n = 3. It 
is true for n — 2 and whenever A has positive matrix off-diagonal elements in an 
eigenbasis of B [BMV]. Moreover, it is true for sufficiently small A [FW]. The 
measure is known to have support in the convex hull of the eigenvalues of B and 
to have a positive atomic part, supported by the eigenvalues themselves. A review 
of the existing partial results, together with some new ones, is in preparation [FW]. 
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Figure Captions 



Figure 1 (for Example 1): state space of C({1,2,3,4}) with £^ah 
embedded. The triangle is convf^^® = 



Figure 2 (for Example 3): section of state space, with £^ah niarked 
(/? = 3, a = 0.3, 7 = 0.05). The point on the vertical axis is uj^^^. 
It is not contained in the convex hull of £ah- 



Figure 3 (for Example 5): state space of M.2 in logarithmic scale; 
inner circle is In^AH, extension by curves is ln£R,(A42), angles 
determine Iu^ee- 



Figure 4 (for Example 5): upper left: the cones Ci, (72,6*3 from 
(5.7). Remaining panels: Iu^'ah + C*!, i = 1, 2, 3. The intersection 
of these sets is shown in Figure 3. 



Figure 5 (for Example 6): left panel: state space of C({1,2,3}); 
right panel: same figure in logarithmic coordinates (same orien- 
tation). Thick line is £^ah, lens shape is £^r., and dashed corners 
together with Ear are the boundary of £^ee- 
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